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We introduce a family of strongly-correlated spin wave functions on arbitrary spin-1/2 and spin-1 
lattices in one and two dimensions. These states are lattice analogues of Moore-Read states of 
particles at filling fraction 1/q, which are non-Abelian Fractional Quantum Hall states in 2D. One 
parameter enables us to perform an interpolation between the continuum limit, where the states 
become continuum Moore-Read states of bosons (odd q) and fermions (even q), and the lattice 
limit. We show numerical evidence that the topological entanglement entropy stays the same along 
the interpolation for some of the states we introduce in 2D, which suggests that the topological 
properties of the lattice states are the same as in the continuum, while the ID states are critical 
states. We then derive exact parent Hamiltonians for these states on lattices of arbitrary size. By 
deforming these parent Hamiltonians, we construct local Hamiltonians that stabilize some of the 
states we introduce in ID and in 2D. 
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I. INTRODUCTION 

The Fractional Quantum Hall (FQH) effect is one 
of the most fascinating phenomena in strongly corre¬ 
lated electronic systems, in which the electrons of a two- 
dimensional electron gas subject to a strong magnetic 
field form an incompressible quantum liquid supporting 
fractionally charged quasiparticle excitations. The un¬ 
derstanding of this paradigm of topological order was in 
large part made possible by the discovery of analytical 
wave functions, such as the Laughlin’s wave function 
describing the electrons in a partially filled Landau level. 

Since its discovery in 1987 [5], one FQH state has at¬ 
tracted a lot of attention : Unlike the states at filling 
factors with odd denominators, the ;/ = 5/2 FQH state 
with electrons occupying the second Landau level cannot 
be explained by a hierarchical construction based on the 
Laughlin’s states [3]. This opens the door to the possi¬ 
bility of electron pairing and emergence of non-Abelian 
quasiparticle excitations. Indeed the leading candidate 
for the description of the = hj^ FQH state is the 
Moore-Read “Pfaffian” state at filling 1/2 [IHS], describ¬ 
ing the wave function of the electrons in the second Lan¬ 
dau level. Moore-Read states have a wave function de¬ 
fined by |1] 

-ipiwi,. . .,Wm) OC n(wi - wj)«pf 

i<j 

( 1 ) 




where Wi are the positions of the particles on the complex 
plane, 1/g is the filling factor and the magnetic length 
has been set to one. They support fractionally charged 
non-Abelian anyons possessing Majorana fermion states 
at zero energy [HIS]. These non-Abelian anyons have 
attracted a lot of attention due to their applications to 
topological quantum computation liiin]. 


There has been a lot of recent interest in finding FQH 
physics in other systems. Indeed, it is a fundamental 
problem to understand the origin and properties of states 
supporting non-Abelian anyons. Finding different sys¬ 
tems exhibiting the same physics is an important step to¬ 
wards an explanation of these phenomena, for which our 
understanding is still far from complete. Moreover, ex¬ 
perimental realizations of FQH states and manipulation 
of their quasiparticles are a challenge. As such it is in¬ 
teresting to search for new possibilities for realising FQH 
physics and variations thereof experimentally. In partic¬ 
ular, lattice models without Landau levels may pave the 
way towards experimental realization and investigation 
of FQH-like states, for example in optical lattices. 

Two main approaches have been followed to recreate 
the FQH effect in lattices. In the first approach, one 
tries to mimic electrons in a magnetic field by replacing 
the fractionally filled Landau level by a nearly flat frac¬ 
tionally filled Chern band and adding local interactions 
[IIHII]. Non-Abelian FQH states where found in such 
lattice models with topological flat bands [IMT^ . 

In the second approach, instead of trying to reproduce 
the interactions between the electrons and a magnetic 
field, the focus is on reproducing the FQH wave func¬ 
tions on lattices. This approach first started with the 
introduction by Kalmeyer and Laughlin of the bosonic 
Laughlin state at filling fraction 1/2 on a square lattice 
[20) . for which a parent Hamiltonian was later derived 
on the torus USUI] and for more general lattices in the 
thermodynamic limit [33H35] (see also [351 HZ] for related 
results). Hamiltonians were also derived for non-Abelian 
chiral spin liquids with excitations with SU{2)k statistics 
[24] ■ but only in the thermodynamic limit. For k = 2 this 
corresponds to a bosonic lattice Moore-Read state at fill¬ 
ing fraction 1 on a spin-1 lattice [28j. Moore-Read states 
of bosons have also been considered on one dimensional 
lattices, were parent Hamiltonians have been obtained 
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P5] . Other filling fractions and Moore-Read states on 
spin-1/2 lattices have not been introduced before. 

In the continuum, a useful description of FQH states [1] 
uses wave functions expressed in terms of correlators of 
the related edge Conformal Field Theory (CFT). This de¬ 
scription was extended to lattice systems in one and two 
dimensions [ 11130 ] by writing the wave functions of spin 
systems as CFT correlators. The states obtained can be 
seen as an infinite dimensional version of Matrix Product 
States (MPS) which gives a unified treatment of ID and 
2D lattices systems. An alternative way of implementing 
MPS ideas to 2D FQH systems has been developed in 
m- For the Laughlin states, this construction of wave 
functions written as CFT correlators provides states that 
are close to, but not exactly the same as the Kalmeyer- 
Laughlin states on a lattice of finite size, but that become 
the same in the thermodynamic limit |23j . This modifi¬ 
cation of the lattice construction has made it possible to 
construct exact parent Hamiltonians for strongly inter¬ 
acting lattice spin systems of arbitrary sizes [H [5M55] , 
including topological FQH states such as Laughlin states 
of hardcore bosons and fermions [25] . It was shown in 
some cases that these states could be stabilized by a local 
Hamiltonian [31 [31 [31 [31 ■ This description has been 
applied to find a Hamiltonian for the SU(2)2 spin mod¬ 
els but the focus was on the one dimensional spin chain 
[32| and the corresponding two dimensional non-Abelian 
FQH state has not been analyzed so far. 

In this paper we fill this gap by extending the con¬ 
struction of lattice wave functions from CFT correlators 
to non-Abelian FQH states. Using this approach we con¬ 
struct a family of lattice versions of the Moore-Read state 
at filling fraction 1/q. This family of states allows us to 
interpolate between the continuum limit, where all states 
become continuum Moore-Read wave functions, and a 
lattice limit. The states are defined on arbitrary spin- 
1/2 and spin-1 lattices in one or two dimensions, where 
the value of each spin can be mapped to an occupation 
number of bosonic (odd q) or fermionic (even q) particles. 
We provide numerical evidence that the states are criti¬ 
cal states in one dimension and that in two dimensions 
they are topologically ordered states. It is shown that 
states defined on a spin-1/2 square lattice (for q = 2) and 
states defined on a spin-1 square lattice (for q — 1) have a 
topological entanglement entropy which is constant along 
the interpolation to the continuum. This suggests that 
the lattice states have the same topological properties as 
Moore-Read states in the continuum. 

Since we use an approach based on the wave functions, 
it is then relevant to ask whether these states are ground 
states of physical Hamiltonians. Using properties from 
CFT, we derive parent Hamiltonians for the wave func¬ 
tions in the lattice limit. These parent Hamiltonians have 
long-range interactions, are exact on any lattice of arbi¬ 
trary size and we find numerically that they have a non¬ 
degenerate ground space for g < 2. By deforming these 
parent Hamiltonians, we then show numerical evidence 
that the state at filling fraction 1 on a spin-1 square lat¬ 


tice can be stabilized by a local Hamiltonian in one and 
two dimensions, while the state at filling fraction 1/2 on a 
spin-1/2 square lattice can be stabilized by a local Hamil¬ 
tonian in one dimension, which is a first step towards an 
experimental realization of these states. 

The paper is organized as follows : In Sec. |H| lattice 
Moore-Read states are defined from correlators of con¬ 
formal fields. It is shown in Sec. m that these states re¬ 
duce to Moore-Read states of particles in the continuum. 
Properties of the states are computed in Sec. IV where 
evidence that the states are critical in one dimension and 
that the topological properties remain the same along 
the interpolation between the continuum and the lattice 
states in two dimensions is presented. Parent Hamilto¬ 
nians are derived in Sec. |V| Finally Sec. |VI| provides nu¬ 
merical evidence that some of these Hamiltonians can be 
deformed into local Hamiltonians stabilizing the lattice 
states. 


II. LATTICE STATES FROM CORRELATORS 
OF CONFORMAL FIELDS 

Let us consider a lattice with N sites at positions Zj, 
j G {1, 2,..., N} in the complex plane. We will refer to 
one dimensional models when the Zj are restricted to the 
unit circle and two dimensional models otherwise. Let 
a be the average area per site (in ID a is the average 
distance between two sites on the unit circle) and rj = ^ 
a positive real number. The local basis at site j is |nj), 
where rij is an integer that will be interpreted as the 
number of particles at site j. We consider two classes of 
models, labelled by S' S {j, 1} : 

• Models with S = | are defined on a Hilbert space 
of size 2^ : the two states in the local basis are nj G 
{0,1} and correspond to the absence/presence of a 
particle at site j. These states can be expressed in 
terms of spin ^ variables at each site : Si = 2ni — 1. 

• Models with S = 1 are defined on a Hilbert space 
of size 3-^ : the local basis is rij G {0,1, 2}, which 
we interpret as the presence of 0/1/2 particles at 
site j. These states can be expressed in terms of 
spin 1 variables at each site : Si = rii — 1. 

In general, a wave function defined on one of these two 
spaces will have the form 

\ip)s= ^ il){ni,...,nT^)\ni,...,nT^)s. ( 2 ) 

ni ,...,nN 

In the following, we will consider wave functions for which 
the coefficients ... ,niy) can be expressed as the 

correlator of some conformal operators. Let us introduce 
the operators 


( 3 ) 
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where <^(z) is a chiral bosonic field from the c = 1 Confor¬ 
mal Field Theory (CFT), x is a Majorana fermion field, 
: ... : denotes normal ordering, g is a positive integer and 
Srij is 1 if Uj = 1 and 0 otherwise. We also define a phase 
coefficient 


(4) 

where the are phase factors to be specified. 

We propose to consider the wave functions defined by 
the correlator of the previous operators : 

ipsini,... ,nN) oc (V„i(2i)... V„„(zAr)), ( 5 ) 

where 

V^^{Zj)=Ur.,{zj)Vn^{z,). ( 6 ) 

These wave functions are labelled by three parameters : 
q, 77 , S and will be referred to as the {q,r])s CFT states. 





FIG. 1. Illustration of a square lattice on the complex plane 
in the continuum limit (77 —>■ 0 , —>■ 00 ) and in the lattice 
limit(77 —>■ 1). At each site there can be 0 (blue circle), 1 
(blue disk) or 2 (red disks) particles. The interpolation is 
performed by fixing the number of particles M = and by 
varying rj — ^ between 0 and 1, which changes the number 
of lattice sites per particle between infinity and q. 

Let us now evaluate the correlator. Note that we do 
not need to add a background charge in this setting. The 
correlator is zero unless = r/N/q. This condi¬ 

tion fixes the total number of particles in the system to 
M = ''7/9 is therefore the lattice filling 

fraction and 77 is a parameter which interpolates between 
the continuum limit (77 —> 0, > 00 , number of particles 

conserved) with infinitely many lattice sites per particle 
and the lattice limit (77 = 1) at which the lattice filling 
fraction is equal to 1 /q, which corresponds to the Landau 
level filling fraction in the FQH effect (Fig. [^. Let us 
write = 1 if ni = M and = 0 otherwise. The 
evaluation of the correlator yields |33 


7/is(77i, . . . ,7ZAr) OC^ ^[(Zj - 
i<j 

X PL ' ^ 






(7) 


where fwizi) = ^iY\j(^i){zi — zj) ^ and the Pfaffian is 
evaluated at the coordinates where tt.^ = 1. The Pfaf¬ 
fian is antisymmetric, so these are bosonic states when 
q is odd and fermionic states when q is even. Note that 
this formula defines different states for S = 1/2 than for 
5 = 1, the difference being that in the second case the 
rii can take the value 2. The states {q, 77 ) 1/2 are therefore 
projections of the states {q, 77)1 onto the Hilbert space al¬ 
lowing only for single occupancy at each site, while states 
with 5=1 and q odd (resp. even) can have sites with 
two bosons (resp. fermions of different types). The state 
(g = 1,77 = 1 ) 1/2 is trivial since the number of particles 
is fixed to N and there are N sites, so in the following 
we restrict ourselves to states (9,77)1/2 with q > 2 and 
(9,77)1 with 9 > 1 (Table |l]). 


TABLE I. First ( 9 , 77)3 CFT states 

Lattice limit 77 = 1 

S 

q 

1 

2 

1 

1 

X 

bosonic 517(2)2 (9 = 1,77 = l)i 

2 

fermionic (9 = 2,77 = 1 ) 1/2 

fermionic (9 = 2,77 = l)i 


Continuum limit 77 —>■ 0 ^ 


5 

q 

1 

2 

1 

1 

X 

bosonic Moore-Read 

2 

fermionic Moore-Read 

fermionic Moore-Read 


III. THE CFT STATES BECOME 
MOORE-READ STATES IN THE CONTINUUM 
LIMIT 

In this section we consider a two dimensional lattice 
defined on a disk V of radius i? —>■ 00 and show that the 
CFT states we have introduced reduce to Moore-Read 
states of particles in the continuum, that is Eq. Q , when 
77 —0, TV — 00 and the number of particles M is fixed. 
We restrict ourselves to lattices where the area per site 
Qi is constant equal to a, but the derivation remains true 
for any lattice if we make 77 position dependent [ 25 j . 

Let us first compute Yii fN{zi)^‘ ■ Notice that I/at] = 
exp(-I]/(^ 0 ^ 1 n(|z/ - Zj\)) and since 77 = ^, in the 
continuum limit this sum can be replaced by an integral 
Jj,lu{\zi — z\))dz'^/2 tt. This integral evaluates, in the 
thermodynamic limit, to + constant [1], so that 

fN^oo{zi) OC 

where gi = Im — Zj)^ is a real number. It 

was found numerically in Ref. [2S] and Ref. [53] that this 
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formula was an accurate approximation even for moder¬ 
ately large N. We thus get that 

Y[fN^oo{zir oc (9) 

In the rest of this section we set the phase factors such 
that = e*®' to get rid of the overall gauge factor, which 
does however not change properties like the particle- 
particle correlation function and entanglement entropy 
of the state. 


A. Continuum limit of the S = ^ states 

Let us now write the complete wave function in the 
continuum limit : 


Since we interpret the state |2) as the presence of two 
particles, positions Zi where rii = 2 are listed with two 
different indices r in the set {wr}. We now write the 
wave function in the basis given by the sets {rcr-}. 

As a starting point, observe that 

M n 

g-|»dV4 ^ g-n.P.lV4_ (12) 

r=l 2=1 


We will now prove that 


r<s 


sYPi 


Wr — Ws 


i<j 



ips(.ni,...,nN) oc(5„ 

i<j 

X Pf„ ' ^ 


4-2 -C-J 


g-3 Ei 

( 10 ) 


where the gauge factor has been set to one. It is not 
straightforward to take the continuum limit in this basis, 
since one has to define the limit of the Hilbert space on 
which the wave functions are defined. However, since the 
number of particles is conserved, we can rewrite the wave 
function in the basis spanned by the positions ici,..., wm 
of the particles. For ^ = 1/2 there is at most one particle 
per site so the wave function can be simply expressed as 

. . .,Wm) oc - Wj)'^Pi 

i<i 

( 11 ) 


e-iT.iMY 


where the Wi are restricted to positions in the lattice. In 
the limit of infinitely many lattice sites per particle the 
lattice becomes a continuous plane and the positions Wi 
become positions in the plane. This state then coincides 
with the Moore-Read state 0- The number of particles 
on the lattice is M = ri— = t/- —, so if the flux is 
then we can express 1 = This explains that this 

quantity corresponds to the filling fraction in the con¬ 
tinuum, defined as the number of particles per magnetic 
flux. 


B. Continuum limit of the S' = 1, q = 1 state 


For S = 1, g = 1, the state also has the form (10), how¬ 


ever since the rii can take the value 2, it is not straight¬ 
forward to take the continuum limit. We first have to 
write the wave function in the basis spanned by the po¬ 
sition of the particles. For a basis element |ni,..., tiat), 
let Wr, r G {1,... ,M} be the positions of the particles. 


Suppose first that only the first site is doubly occupied. 
Then, 


- Ws) = (zi - Zi) 0 (zi - 0 {Zi - Zj)'^^' 

2 <i<j 

(14) 


r<s 


1<2 


Moreover, 


Pf 


Wr — Ws 


1 


Zi - Zi 
1 

Zl - Zi 


1 


-Pf„ 


Wr — Wg 

1 


(15) 

(16) 


where we have used in the last line that all other positions 
have at most one particle. Combining (14) and (16) we 
see that the (zi — z\) term cancels, so that (13) holds. 
The same derivation with a recursion on the number of 
sites with two particles proves that (13) holds for all basis 
elements. 

Putting together (13) and (12), the wave function can 

1 


therefore be written as 


ipl ^(wi,..., Wm) oc 0(u;i - Wj/^Pf 

i<j 


Wi — Wj 


(17) 


where the Wr can be repeated twice to allow for states 
with double occupation, in which case this expression 
does not vanish because of a cancellation between the 
Jastrow factor and the Pfaffian. If we now take the con¬ 
tinuum limit, the positions of the particles can be any¬ 
where on the plane and this becomes the bosonic Moore- 
Read state ([^ at filling fraction 1, which also does not 
vanish when two particles are at the same site. 

One may ask what happened to the doubly occupied 
site. In the continuum the ensemble of states with two 
particles at the same positions has measure zero com¬ 
pared to states with at most one particle at each posi¬ 
tion, therefore they are irrelevant and do not contribute 
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to the wave function. Similarly, when we take the contin¬ 
uum limit of the lattice, states with two particles at the 
same site do not contribute to the wave function, as we 
will now explain in more details. Intuitively, this comes 
from the fact that there is a finite number rii = M 
of particles, so configurations with at least one doubly 
occupied site are rare when the number of sites goes to 
infinity. More quantitatively, let us denote Pm the num¬ 
ber of basis elements |ni,..., un) satisfying ^ = M 

with m doubly occupied sites. Pq = (^) is the num¬ 
ber of basis elements with no doubly occupied sites and 
Pm = m\(M-2m^{N+m-M)r the continuum limit, M 
and m are fixed and N goes to infinity, so Pm/Po 
K/N™' where K is a constant and Y^m=iPmlPo 0. 
This shows that the number of basis elements with at 
least one doubly occupied site is small compared to the 
number of basis elements with no doubly occupied sites. 
Let us now observe that the wave function has to be nor¬ 
malized by a factor Q = ■ ■ ■ ,njv)p. We 

can decompose Q between the basis of states with zero 
doubly occupied sites and the basis of states with at least 
one doubly occupied site : 


C. Continuum limit of the S = 1, q>2 state 


When q>2 the state can be written as 

... ,nAr) cx ]^(zi - ... ,njv) 

i<j 

( 21 ) 

We have already derived the continuum limit of the state 
on the right at g = 1 and the remaining factor can be 
expressed as 


i<j r<s 

where 


g{Wr,Ws) 


Wr — Ws if Wr ^ Wg 

1 otherwise 


(23) 


In the continuum limit, the wave function can therefore 
be written as 


Q= Y1 


(18) 

ni , 



= E 

5j2ni=M\i^{ni,... ,nAr)p, 

(19) 

ni , 



II 

M 


..,n^)|2, 

E rzi^M 

ni,...,n7v 

E rii^M 



( 20 ) 


where the first sum contains Pq elements and the second 
sum contains Y^m=i Pm ^ elements. We now ob¬ 
serve that elements appearing in these two sums are of 
the same order. Indeed, let us take a configuration with 
M particles at positions Wr such that 2 particles are at 
the same site (wi = W 2 ), corresponding to an element 
ip^ni,... ,njv) in the second sum. Then by slightly mov¬ 
ing one of the particles {wi = Wi + S) we obtain a new 
configuration that is in the first sum, where all particles 
are at distinct positions. Since the number of particles 
is fixed, continuity of the coefficients of the wave func¬ 
tion 0 implies that the new element obtained is close 
to ..., un)- Therefore elements in the second sum 
are of the same order as elements in the first sum and 
since there are X)m=i Pm ^ ^0 elements in the second 
sum, in the continuum limit this sum does not contribute 
to the normalization factor Q. This shows that the con¬ 
tribution of configurations with two particles at one or 
more sites can be neglected in the continuum limit. 

Note that the previous derivation also shows that the 
(g = l,r] = 1)1 CFT state in the thermodynamic limit 
is equivalent to the spin 1 non-Abelian chiral spin liquid 
introduced in Ref. [55], but the two states are different 
on finite lattices. 


. . . , Wm) OC g{Wr - Wg)'^ 
r<s 


X Pf 



J|e-|-dV4^ (24) 


When no two particles are at the same site, this wave 
function is the same as the (g, g)i /2 CFT state in the con¬ 
tinuum limit. As explained in Sec. |IIIB[ configurations 
with two particles at the same site do not contribute to 
the wave function in the continuum limit, therefore the 
{q,g)i CFT states for g > 2 have the same continuum 
limit as the (g, 77 ) 1/2 CFT states, which is the Moore- 
Read state at filling fraction 1/q. 


D. One dimensional continuum limit 

So far we have focused on two dimensional states. In 
the one dimensional setting, when Zj = , the same 

results enable us to perform an interpolation between the 
lattice and the continuum. The only difference with the 
2D case is that it is now possible to compute analytically 

fN{zi) OC (25) 

so the wave functions in the continuum can be expressed 
as 


ips{wi, ■. . ,U>m) OC ^(■11;, - 
Kj 



(26) 


which is a one-dimensional version of the Moore-Read 
state. 
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IV. PROPERTIES OF THE GET STATES 

In two dimensions, the Moore-Read states in the con¬ 
tinuum are topological states which support non-abelian 
quasi-particle excitations. It is of high interest to check 
whether the lattice CFT states we have introduced share 
these properties. In one dimension, we expect that the 
CFT states display critical behaviour related to the con¬ 
formal operators used to construct the wave function. In 
this section we focus on the states with q—1 and 9 = 2 , 
and numerically compute some of the properties of the 
states we have introduced. 


A. One dimensional critical states 

We now look at one-dimensional chains such that Zj = 
^ 2 -n:ij/N ^ Since we will find local Hamiltonians for the 
states (9 = 1 , 77)1 and (9 = 2 , 77 ) 1/2 in one dimension in 
the lattice limit in Sec.|VI[ we focus on these states. First 
we compute the Renyi entropy = — InTr where 
Pi is the density matrix of the CFT states restricted to a 
subsystem of size L of the chain. This computation can 
be performed by using a Metropolis-Hastings algorithm 
with two independent spin chains [301 13HH3D] . 

The results are shown in Fig. The entropy scales 
logarithmically with the size of the subchain for all values 
of 77 . Moreover, the scaling is approximately the same for 
different values of 77 and fits of the form 

Sl'’= + b (27) 

yield a value of the central charge c approximately equal 
to 1.36 for the (9 = 1 , 77)1 CFT states (the main source 
of errors is here the finite size of the lattice considered). 
This value is in agreement with the value of 1.395 found 
for the state in the lattice limit in Ref. [32], where it was 
also shown that a value of 1.5 for the central charge, 
as expected for the SU{2)2 WZW model, could not be 
excluded. For the (9 = 2 , 77 ) 1/2 CFT states we find a 
value of 0.98 which is compatible with a central charge 
equal to 1. 

Another quantity that can be computed using Monte- 
Carlo techniques is the particle-particle correlation func¬ 
tion Cl = (uittl) — {ni){nL). Results in Fig. confirm 
that the states are critical since the correlation func¬ 
tions decay polynomially with the distance L. For the 
(9 = 1,77 = 1)1 CFT state, the critical exponent is found 
to be 0.70, which is in agreement with the value of 0.69 
found in Ref. [33] , where it was observed that such a value 
can be influenced by a multiplicative logarithmic correc¬ 
tion which could explain the difference with the expected 
value of 0.75 [3TH33] . Moreover for the (9 = 2 , 77 ) 1/2 
states at different values of 77 , the correlations are very 
close once rescaled by a factor of 1 / 77 ^, which confirms 
that properties of the state do not change along the inter¬ 
polation. For a Tomonaga-Luttinger liquid, the expected 



ln{sin(7rL/A'’)) 

(a) 



ln(sin(xL/A^)) 

(b) 


(2) 

FIG. 2. Renyi entropy ' of a subsystem of L consecutive 
sites for the ID (9 = 1 , 77)1 CFT states (a) and (9 = 2 , 77 ) 1/2 
CFT states (b) for different values of 77 . The number of 
particles M = vN/q is fixed so the sizes of the chain are 
N = 40,80,160,320 for 77 = 1,1/2,1/4,1/8 respectively. The 
lines are linear fits of the points for 77 = 1/8 (blue) and 77 = 1 
(red). 


behaviour of the particle-particle correlation function is 

mi 


Cl = 


Acos{2LkF) 


K 


|sin(5|^)^ 


\2K 


2^2 


(28) 


where K is the Luttinger parameter, kp = 7777/9 is the 
Fermi momentum and A is a non-universal constant. For 
the (9 = 2 , 77 ) 1/2 state we find a good agreement of this 
formula for K = 0.494, A — 0.123. This suggest that this 
state in one dimension is well described by a Tomonaga- 
Luttinger liquid with central charge c = 1 and Luttinger 
parameter K = 0.5, which corresponds to the properties 
of a free-boson CFT with radius as was the case for 
the corresponding one-dimensional Laughlin state [35]. 
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- 0.04 

- 0.05 


sin{7rL/N) 

(b) 


FIG. 3. Rescaled correlation function 'q~^CL as a function 
of the distance between the sites for the ID (g = CFT 

states (a) and (q = 2, r;)i /2 CFT states (b) for different values 
of 77. The number of particles M = rjN/q is fixed so the sizes 
of the chain are N = 40,80,160,320 for rj = 1,1/2,1/4,1/8 
respectively. The data for 77 = 1 is shown in the insets in log- 
log scale, conhrming the polynomial decay of correlations, and 
the line in the insets is a linear fit yielding critical exponents 
0.70 {{q = 1 ,77 = 1)1 state) and 1.02 {{q = 2 ,77 = 1 ) 1/2 state). 
In (b), the line is a ht of the form Eq. (28 1 with parameters 
K=0.494 and A=0.123. 


B. Two dimensional topological states 

In the continuum, the Moore-Read state at filling frac¬ 
tion 1/q has a topological entanglement entropy of [IS] 

70 ( 9 ) = ^ln(4g). (29) 

To compute the topological entanglement entropy 
(TEE) of the CFT states, we map a. x Ly square 
lattice on the cylinder to the complex plane by choos¬ 
ing the positions of the lattice sites to be Zj = 

^2-7r{{xj-L^-l/2)+iyj)/Ly ^ -^^liere Xj G . ,Lx}, yj G 

{l,...,Ly} (Fig. 1^. We then cut the cylinder in two 
halves and compute the Renyi entropy of the first half. 
The size along the cut is Ly and we use the behaviour of 


lm{zj) 


Re(2;/) 




-*■ 




FIG. 4. The mapping from a lattice on the complex plane 
to a cylinder. To compute the topological entanglement en¬ 
tropy of the state, the cylinder is cut into two halves and the 
Renyi entropy of the first half is computed using a Metropolis- 
Eastings algorithm. The topological entanglement entropy is 
then extracted by varying the size Ly of the cylinder (301. 


the entanglement entropy [IQ Wf\ 


s''^l=aLy-^ 


(30) 


to extract the topological entanglement entropy 7 (Fig. 
[^. The results for the states {q = 1 , 77)1 (resp. {q = 
2 , 77 ) 1 / 2 ) are in agreement with the topological entangle¬ 
ment entropy of a Moore-Read state at filling fraction 
1 (resp. 1/2), 7 o(l) = 5 ln(4) « 0.69 (resp. 70 ( 2 ) = 
4 ln( 8 ) « 1.04). Moreover the value of the TEE does 
not change with 77 , so topological properties of the states 
remain the same along the interpolation between the con¬ 
tinuum and the lattice limit. 

We observe however that the state {q = 2 , 77)1 has a 
TEE close to zero and different than 70 (2). The TEE 
does not stay constant when 77 is changed, which is com¬ 
patible with the expectation that its value is 70 ( 2 ) in the 
continuum limit : there must be a phase transition along 
the interpolation between the continuum and the lattice. 
The states with S' = 1, g > 2 can therefore define distinct 
lattice states from the states with S = 1/2, q > 2, while 
having the same continuum limit. 

















































Ly 



Ly 


FIG. 5. Linear behaviour of the Renyi entropy with the size 
of the cut Ly for the (g = 1, g)i (a), (g = 2, ??)i /2 (b) and (g = 
2,?7)i (c) CFT states on a Lj; x Ly lattice. The topological 
entanglement entropy of the Moore-Read states at filling 1 
(a), 7o(l) ~ 0.69, and at Hlling 1/2 (b,c), 7o(2) Rs 1.04, are 
indicated with a red arrow. The values of g are 1,1/4 and 1/8 
and the corresponding sizes L^ are respectively 12, 16 and 16. 
The insets are enlarged views confirming that the topological 
entanglement entropy stays the same when g is varied and 
that its value corresponds to 7o(l) (resp. 7 o(l/ 2 )) in the first 
two cases, while the topological entanglement entropy of the 
(g = 2,g)i CFT state is close to zero in the lattice limit and 
close to 7 o (1/2) in the continuum limit. 


V. PARENT HAMILTONIANS 

So far we have considered wave functions for lattice 
states. It is also relevant to ask whether these states are 
ground states of some Hamiltonians and whether these 
Hamiltonians can be realized in nature or implemented in 
experiments. We turn now to the construction of parent 


Hamiltonians for which the CFT states are ground states. 
The CFT states are constructed from correlators of con¬ 
formal fields. As has been shown in Ref. [35], this enables 
the construction of parent Hamiltonians from null helds 
of the considered CFT. Null fields are fields such that 
when inserted in a correlator of primary helds, the eval¬ 
uation of the correlator gives zero. The procedure is as 
follows : 

1. Find null helds labelled by a and acting at 

position Zi of the considered CFT. 

2. The vacuum expectation value of a product of pri¬ 
mary chiral conformal helds is zero if one of the 
helds is a null vector, therefore we get equations 

(V„,(zi)...x“(zi)---V„„(zw)) = 0. (31) 

3. These equations are rewritten in the form A“|?/;) = 
0, where A“ are operators acting on the degrees of 
freedom of the wave function. 

4. ■ A“1A“ is then a positive semi-dehnite Hermi- 
tian operator annihilating the wave function. 

In this section we apply this procedure to construct par¬ 
ent Hamiltonians for the CFT states in the lattice limit. 
In the rest of this work the phase factors are hxed to 
^1 = 1 - 


A. Parent Hamiltonians for the SU{2)2 
{q — l,r] — l)i CFT state 

The special case of the (g = l,g = l)i state has a 
wave function constructed from the spin 1 primary helds 
of the SU{2)2 Wess-Zumino-Witten conformal held the¬ 
ory (35] US] and has already been considered partially in 
Ref. [35], where, however, the focus was on ID systems. 
This SU{2)2 symmetry can be used to construct parent 
Hamiltonians invariant under SU{2) transformations. 

As shown in Ref. [35], the null helds in this case can 
be parametrized by 

X“(z.)= (32) 

where repeated indices are summed over, (/?„ is a chiral 
spin 1 primary held, J^i are the —1 modes of the SU{ 2)2 
current operators and 

= Isab - (33) 

where eabc is the Levi-Civita symbol and are the spin 1 
operators acting on site i. These operators can be written 









in the spin basis at site i as 


9 


.X _ 

* 72 


^0 1 0 ^ 

fv _ A 

^0 -i o\ 

1 0 1 

i 0 —i 

Vo 1 q) 

V 2 

yO i 0 J 


1 0 0 
tf = h I 0 0 0 
,0 0 - 1 , 


(34) 


We then exploit that the correlator with a null field 
inserted is zero : 


{ifm (.Zi)... X^'iZi) ■ ■ ■ {zn)) = 0. 


(35) 


The Ward identity enables to transfer the action of a 
current operator to fields at other positions : 

■ ■ ■ {.JtiiPni){Zi) . ..ipri„{.ZN)) = 


N 


t? 


y" - - - {iPnAzi)---^ni{Zi)---^m,{zN))- (36) 

• * 7 '.-\ — r • 




Zi - Zj 


Using (36 1 it is possible to rewrite (35) as 
A“|^)=0, 

where 

N 


K=T, (^'“’)A«t‘. 


(37) 


(38) 


N 


= E 


UJ. 






(39) 


B. Parent Hamiltonians for the {q,r] = l)i, q>2 
CFT states 


For q>2, the {q, rj = l)i CFT states do not display an 
SU{2) symmetry. However it is still possible to find null 
fields and construct operators annihilating the wave func¬ 
tion. Let us define the operators G^{z) = xiz)e^’‘'^‘^G) ^ 
J{z) = -^d(j){z). We use the following q + 1 null fields : 


X^iw) 

p=0,l,...,(ir-2 


dz 


1 


27ri (z — w)P 


G+{z)V^{w), 


(41) 


x-y-) = 


1 


(z — w)l 

1 


-G+(z)V^{w) 


(z — w) 


X«(w) = 


1 


dz 1 

27ri z — w \_{z — w)i 
^.—^qJ{z)V2{w), 

ZTTZ Z — W 


V2{w) , (42) 

-G+{z)V,iw) 

(43) 


where Eq. (41) gives 9 — 1 null fields x^{'^) with p running 
from 0 to q — 2. Let us define the operators d\ d to be 
creation and annihilation operators (bosonic for q odd, 
fermionic for q even) acting between states | 0 ) and | 1 ), 
and d'^ , d' to be creation and annihilation operators act¬ 
ing between states |1) and |2). The number of particles 
at site i is thus Ui = -|- 2n[^\ where = d'^d and 

= d'^'d'. Using the previous null fields, we derive the 
following operators annihilating the wave function (see 
Appendix [A] for the complete derivation). 


where we have defined uii^ = This can be used to 

J Zi—Zj 

construct a parent Hamiltonian H = A)*^A“, which 

gives 

N \ ^ ( ^ 

= 2 E A g E I 2 E 

\ k{^i,i) 


i 

A? = U 

Af“^ = W +rif\ 

* ^ [zi- zU 9 -i * * 

i(#i) ^ 

A? = E . ^ . d,d!^ _ y 


(44) 

(45) 

(46) 


i_ d d!^ - V 

(z. - z M ^ ^ 2^ z - z ^ ^ ■ 

3{¥^^) 


1 _ /1 -l \ 

“ 6 E + E ( ^^*k^i3 - 2^rk^z] ) 

i^j i¥=3¥=k ^ 


(47) 


i+i+k 


This leads to a three-body Hamiltonian annihilating 
the wave function 


This Hamiltonian is SU(2) invariant and numerical di- 
agonalization on small systems confirm that it has the 
(g = 1, 77 = 1)1 CFT state as a unique ground state. This 
Hamiltonian is similar to the one obtained in Ref. [53] 
for the spin 1 non-Abelian chiral spin liquid state intro¬ 
duced in Ref. [28]. However the Hamiltonian we just 
constructed is valid for any choice of lattice and not only 
in the thermodynamic limit as is the case in Ref. 


N q 



i—\ a =0 



(48) 


where the last term fixes the number of particles. 

This Hamiltonian annihilates the wave function (g, 77 = 
1 ) 1 , however we find numerically that the ground space 
of this Hamiltonian is degenerate when g > 3 and that 
the degeneracy does not depend on the number of sites. 
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Other simple null fields of the theory constructed with 
the same current operators do not lead to operators act¬ 
ing on the wave function that would reduce this degen¬ 
eracy. 


C. Parent Hamiltonians for the {q,ri — l)i/ 2 ! > 2 

CFT states 


We can now use the previous results to construct par¬ 
ent Hamiltonians for the (g, 77 = 1 ) 1 / 2 , <? > 2 CFT 
states. They are projections of the {q,r] = l)i CFT 
states in the subspace allowing only for single occupa¬ 
tion at each site, that we denote Let us also define 
the Hilbert space spanned by basis elements contain¬ 
ing at least one site with two particles. We will now 
project the operators annihilating the {q,r] = l)i CFT 
states onto in order to get operators annihilating the 
{q,r] = 1 ) 1/2 CFT states (see Appendix |B for the de¬ 
tailed derivation). We start by multiplying xne operators 
Af obtained previously on the left by d' : these operators 
continue to annihilate the (< 7,77 = l)i CFT states. Since 
d'df = we get operators A°, d'A),..., d'A®“^ that 
act on and are zero on so they also annihilate 
the (( 7,77 = 1 ) 1/2 CFT states. We can then use the fact 

that d'Af^ = + d'i annihilates 

the (g ,77 = l)i wave function to replace the operator 
d' in d'A« by -Y.h{^i) The resulting 

operator d'A^ then acts separately on 'H} and We 
finally get operators annihilating the ( 5,77 = 1 ) 1/2 wave 
function : 


VI. LOCAL HAMILTONIANS 

The parent Hamiltonians we have derived involve three 
or five-body interactions between all sites on the lat¬ 
tice. These Hamiltonians would therefore be very dif¬ 
ficult to implement in experiments. However in some 
cases it has turned out that states constructed from cor¬ 
relators of conformal fields had very high overlaps with 
ground states of local Hamiltonians [55113011231 [3^ . This 
has lead to a protocol to implement one of these states 
in experiments [26l [49]. In this section we show that 
there is a local Hamiltonian for which the ground state 
is close to the {q = 1,77 = l)i CFT state in one and in 
two dimensions and that this result is also true for the 
(g = 2,77 = 1 ) 1/2 CFT state in one dimension. 


A. Local Hamiltonians for the (g = 1,77 = l)i CFT 
state 


In one dimension the case of the (g = 1,77 = l)i CFT 
state was studied in Ref. [35]. It was shown that this 
state has a high overlap with the ground state of the 
bilinear-biquadratic spin 1 Hamiltonian 


rr(l) _ 
^ID — 


N 

[cos (/3)t“t“+i + sin (^) > (53) 


with periodic boundary conditions, when /3 = 0.3213. 
Note that this Hamiltonian includes the 2-body terms 
present in the parent Hamiltonian (40). 


TABLE II. Terms in the Hamiltonian H 2 D and coefficients 
obtained after nnmerical optimization on a 4 x 4 lattice. 


p=l,...,q-2 


i 

(49) 

Operator 


Configuration 

Coefficient 


(50) 


0 

a 0 0 0 

0 ) 1 

b) 0.6227 

2^ iz- - z ^ ’ 

1 J „(1) 

j.a j.a 

Li Lj 

0 

9 0 9 0 

(z- - z )g^ ’■ 

+ y y 7 \ 


{ttY? 

0 

0 

0 0 

c) -0.1762 

d) 0.3226 

j{^i) h{^i) 

z,-z^ 

(51) 

laiaibib 

Cj Lj tj L]^ 

0 

9k 0 0 0 

I " 

A a a 0 

i i k 

e) 0.46377 

/) 0.0208 


This leads to a five-body Hamiltonian 



As in the previous case, this parent Hamiltonian has 
a single ground state only when g = 2 and this ground 
state is the (g = 2,77 = 1 ) 1/2 CFT state. 


We now study the two dimensional case and build a lo¬ 
cal Hamiltonian from the parent Hamiltonian (40). The 
operators A“ contain 2-body interactions between sites 
i and j. We cut these operators by keeping only terms 
for which the sites i and j are nearest-neighbours on the 
square lattice. This leads to a local Hamiltonian with 
three-body interactions. In addition to these terms, we 
include the two-body interactions between next-nearest 
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neighbours present in the parent Hamiltonian. All six 
terms included in our trial Hamiltonian, that we denote 
H 2 Dj are shown in Table [ll| Note that the coefficients 
of these terms in the exact parent Hamiltonian are not 
position-independent. In our local Hamiltonian, how¬ 
ever, we choose them to be position-independent and in¬ 
variant under rotations. 

By exact diagonalization and optimization on these co¬ 
efficients, we find that there is a local Hamiltonian for 
which the overlap |(^/'//|'0cft)| between the ground state 
and the {q = 1,?7 = l)i CFT state on a 4 x 4 square 
lattice on the plane is 97.36%. Considering the size of 
the Hilbert space 3^® « 4 x 10^, this overlap is very high. 
Note that with the same parameters, the overlap is also 
above 98% on a 4 x 3 or on a 4 x 2 lattice. On a cylin¬ 
der geometry, i.e. periodic boundary conditions in one 
direction, the overlap on a 4 x 4 square lattice is 97.21%. 

Compared to the local Hamiltonian found in Ref. [28] , 
which is for a state that is equivalent to the {q = l,r] = 
l)i CFT state in the thermodynamic limit, but different 
on finite lattices, the Hamiltonian we find has less free 
parameters (5 instead of 11) to fine-tune, which might 
make it easier to implement. Moreover, the very good 
scaling with lattice sizes lets us expect that a good agree¬ 
ment will persist on larger lattices. In Fig. [^ we show 
the low-energy spectrum of this local Hamiltonian. This 
figure is compatible with having a gap in the thermo¬ 
dynamic limit, but the limitations on the system sizes 
that we can consider prevent us from making a reliable 
extrapolation. 


a local Hamiltonian with several terms. We find how¬ 
ever numerically that even a smaller number of terms is 
already sufficient to get a good overlap. Specifically, we 
choose to keep only the simplest two-body and the sim¬ 
plest three-body terms, to obtain a local Hamiltonian 
with periodic boundary conditions 

N 

where dl, di are fermionic creation and annihilation 
operators at site i and rn = d\di. For k = 0.274-|-0.052z, 
we find that the overlap between the ground state of 
and the (q = 2, rj = l )i /2 C FT state is 97.71% for 
a chain with 20 spins (Fig. |7(a) I. 



N N 


(a) (b) 



FIG. 6 . Energy difference to the ground state energy for 
the first excited states of the Hamiltonian H 2 D for different 
sizes of lattices. The overlap between the ground state of this 
Hamiltonian and the (g = 1,77 = l)i CFT state is indicated 
below each ground state. 


B. Local Hamiltonian for the ID {q = 2,ri = 1 ) 1/2 
CFT state 

The state (g = 2,77 = 1 ) 1/2 CFT state has a five-body 
parent Hamiltonian. Let us consider the one dimensional 
case. If we cut the A“ operators by keeping only terms 
for which the sites i and j are nearest-neighbours, we get 


FIG. 7. (a) Overlap |{V’CFT|'!/'n)| between the {q = 2,ri = 
1 ) 1/2 CFT state and the ground state of Hamiltonian 
with K fixed to 0.274 -|- 0.052i, as a function of the number 
N of lattice sites. The dotted fine is a linear fit with equa¬ 
tion y = 1.0075 — 0.00152A. If the overlap continues to fol¬ 
low this behavior at larger sizes, it will still be above 85% 
for a spin chain with 100 lattice sites, (b) Overlap per site 
between the same two states. 


In two dimensions, since the SU{2) symmetry is not 
present in this model, cutting the parent Hamiltonian 
leads to a local Hamiltonian with up to five-body inter¬ 
actions with many different coefficients. In addition, the 
fact that the five-body terms involve more sites means 
that each of them stretches over a larger part of the lat¬ 
tice. With the limited lattice sizes that we can consider 
with exact diagonalization, this is problematic because 
the local regions need to be small compared to the total 
size of the lattice (otherwise it would not be expected that 
the same local Hamiltonian would also work for other lat¬ 
tice sizes). This suggests that even if a local Hamiltonian 
that is related to the exact Hamiltonian exists, we may 
not be able to find it with exact diagonalization. Instead 
of cutting the exact parent Hamiltonian, we therefore 
asked whether, by chance, a local Hamiltonian can be 
obtained if we restrict the range of the interactions to all 
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interactions preserving the number of particles on all pos¬ 
sible configurations inside a plaquette of the lattice. Op¬ 
timizing the coefficients in this Hamiltonian, we did not, 
however, find a set of coefficients for these interactions 
for which the ground state of this Hamiltonian is close to 
the {q = 2,r] = l)i /2 CFT state. Whether there exists a 
more complicated, but still local, Hamiltonian stabilizing 
this state therefore remains an interesting open problem. 


VII. CONCLUSION 


We have introduced a three-parameter family {q,ri)s 
of strongly-correlated spin states on arbitrary lattices in 
one and two dimensions. It was shown that these states 
reduce to continuum Moore-Read states of bosons (odd 
q) and fermions (even q) in the continuum limit ry —>■ 0. 
Numerical evidence that these states are critical states in 
one dimension and topological states in two dimensions 
was provided, and the topological entanglement entropy 
was shown to remain the same along the interpolation for 
the {q = and {q = 2,r])ii2 states. Parent Hamilto¬ 
nians of the states in the lattice limit were derived using 
analytical tools from CFT and in some cases it was shown 
that these states could be stabilized by local Hamiltoni¬ 
ans in one and two dimensions. 

There is currently a lot of interest in finding models 
possessing topological properties. Given the complex¬ 
ity of quantum many-body systems, the analysis of phe¬ 
nomena, like e.g. topology, can be greatly facilitated by 
having models in which at least the ground state can be 
found analytically. The results of the present paper show 
that CFT is a valuable tool to derive analytical models, 
also in the context of non-Abelian FQH states, and in ad¬ 
dition can be used as a starting point to identify simpler 
models that are more realistic to realize physically. 
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Appendix A: Operators annihilating the {q,rj = l)i 
CFT wave functions 

We restrict ourselves to the case q > 2, since the con¬ 
struction of the parent Hamiltonian when g = 1 is ex¬ 
plained in |V A[ Note however that the derivation pre¬ 
sented here can also be used to obtain a Hamiltonian 
without the SU{2) symmetry for the (g = 1,?7 = l)i 
CFT state. The CFT states are defined from the opera¬ 
tors Vnj{zj) given in Eq. Q. Additional operators that 


are needed to construct parent Hamiltonians are the op¬ 
erators G±(z) = J{z) = -^d(t){z). The 

operator product expansion (OPE) of G+(z) and J{z) 
with the operators used to build the wave function are, 
for q > 2, 

G+{z)V^^iw) ~ (Al) 

z — w 

J{z)Vr,yw) ~ (A2) 

q z — w 

We first show that a null field, where 


J W 


dz 1 


1 


27rf z — w {z — wy 

qJ{z)V2{w) 


-G+iz)Vyw) 


2TTi z — w 
= - 0 |(w), 


(A3) 

(A4) 


where the contour of the integration is a circle around w 
traversed counter-clockwise. We have that 


= i 

J W 


and 


dz 

1 

1 

2Tri z 

— w 

o 

1 

1 

dz 

1 

1 

2TTi z 

— w 

7 

T 

1 

o 

_ ! 


xe 

+iy/q<j>{z)y(q- 

dz 

1 

'{z-wY-\ 

2Tri z 

— w 

1- 

o 

1 

1 

x. 


xe 

i^<p(z)+i(q-l 

dz 

1 

\fl , 



- + 


-G+{z)Vyw) 


27Ti z — w 


z — w 


xiz)xiw) 


+ (z — w)A{w) 


Xg»V4<^(^)+d9-l)0('i")/v^ 


dz 


1 


-I- A{w) + ... 


2-Ki [ \ (z — wY 


dz 1 
27ri [{z — wY 
dz 1 




ni{w) = (f 

J w 


27rf z 

— w 

dz 

1 

27rf z 

— w 

dz 

1 

27rj z 

— w 

dz 

1 


fv^a^(u;)e*(29-i)'^(“')/vA 


(AS) 


27 Ti z — W - 


[qJ{z)V 2 {w)] (A6) 

V^*a(/)(z)e*(29-i)0(*«)/V9l (A7) 




(A8) 

(A9) 
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which shows that x'^(w) is a null field. Similarly, there 
are more simple null fields 


x’’ i(u>) = I 

J w 


dz 

2Tri 


1 


{z — wY 
1 


^G+{z)VYw) 


{z — w) 


1/2 (w) 


(All) 


Let us now use the fact that by replacing the field at 
site i by a null field, the correlator vanishes : 


0 = ■ ■ ■ X°'iZi) ■ ■ - Vni^izN)). 


(A12) 


We will transform this equation into an equation involv¬ 
ing the wave function by deforming the contour integral 
and moving the operators in the null fields at different 
positions. Let us do it for the null field x'*(w). We will 
use the OPEs as well as the commutation relations : 

V„,.(z,)G+(z) = (-1 )G+i)".-ig+(z)V„,.(z,), (A13) 

where we have used : :: g*/3'/>(“) • = • 

gi/ 30 (^) .. gia 0 (z) . xiz)x{w) = i-YxMxY), which 
adds a minus sign only when rij = 1. 

We then have, starting with the first term involving 

xliw) : 


dz 1 

^i{z-ZiY^^"^ 
dz 1 


(Vni( 2 l) ■ • ■ G+(z)Vi(Zi) . ..Vnt,{zN)), 


= -E 




27ri (z — ZiY 


X (V„i(2:i). ■. G+(z)Vi(Zi) . ■■Vnj,{zN)), 

^ Jz, 27ri (z - z,Y Z - Zj 

X (V„i(zi)... Vi(zj)... Vi(zi)... V„„(zAr)) 
^ /- 27rz (z — ZiY 


j=i+i 


X (V„, (zi )... Fi (z,)... Vi (z, )... V„„ (zjv)), 

Z Zj 

. ^ (_i)(9+i)ElL>f= 

= -(-!)*■' E ^n,.o 

(Zj - ZiY 

X (V„,(zi)... Vi(z,)... Vi(z0 ... V„„(zAr)) 

^ (_l)(.+i)(_l)G+i)EilU.n. 


j=i+l 


(Zj - ZiY 


X (VnYzi)--.Vi(zY... Vl(z,) . . . V„„(zAr)), 

^ (_l)(9+l)EtZ + i"fc 

= - A pyny 


^ , . . . , 1 , . . . , 1 , . . . , 


X Y{q,rj=l)i (ni,... ,1, ■ ■ ■ Y, ■ ■ ■ , n-N), 

^ , ('_1')(9+1) A’fcZ + irifc 


i=i 


(_l)(9+l)AZ)+l"fc 

E 7Z iWii E.o 


i=*+i 


(Zi - Zj)9 


X '0(g',77—; 


(A14) 


where Y(q,r}=i)i is the wave function of the ( 5,77 = l)i 
CFT state. 

Let us now define the creation and annihilation oper¬ 
ators dj,dpd,[,d'i acting on the Hilbert space at site j 
as 

{ 0 rij = 0 

| 0 ) n, = 1 (A15) 

0 rij = 2 


| 1 ) nj=0 

d]\nj) = (-!)(«+!)"7= <j 0 n, = 1 

0 rij = 2 

0 rij =Q 

d'|n,) = (-!)(«+') "7= 0 n, = 1 

. 11 ) ri,=2 

0 = 0 

= n, = l 

0 n-j = 2 


(A16) 


(A17) 


(A18) 


We also define the particle number operators correspond¬ 
ing to these operators as = (Yjdj and = d'^d'j, 


such that rii = + 2n^^\ 

J ' 


3 J’ 


We multiply Eq. (A14) by 


lui,..., rij—I , rij , rijj-i ... 5 rii— 1 , 2 , riij-i . ■ - ? ri ^), (A19) 

and sum over all n^, fc 7 *^ i, to get an expression involving 
the wave function 


E 


1 

{Zi - ZjY 




(A20) 


Let us now look at the second term involving x|(ri;): 
- 9 ■ • ■ J(z)V2{Zi) . . . V„„(zAr)), 


= 'iE 


dz 1 
2Tri z — Zi 


(A21) 

{VnAzi)---Jiz)V2{zY...Vn^iZN)), 

(A22) 
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= E ^'E^(V„,(zi)...F 2 (z,)...V„„(^^)). (A23) 
A/i) 


Let us multiply this expression 

by |ni,... ,ni_i, 2 ,ni+i... ,nAr) = 

Sn' Wi, ■ ■ ■ ,n[,..., uat), which leads after sum¬ 
ming over all uj., k ^ i to 


- E 

., , .V 


(A24) 


Summing the two terms together, we get an operator 
annihilating the wave function : 


A|=E 




1 

{Zi - Zj)' 


M - E 




3(.¥=i) 


Zi — Zi 


(A25) 


The same procedure applied to the other null fields 
gives the following operators annihilating the {q,rj = l)i 
CFT wave functions : 


A° = Erf- 


Af = E ^ 




Af-^ = E 


j{¥=i 


[z,, - ZiY 


did'} + 


_ 1 ^3 i 


(A26) 

(A27) 

(A28) 


Appendix B: Operators annihilating the [q,r] = l)i /2 
CFT wave functions 


To obtain operators annihilating the (g, r] = l)i /2 CFT 
wave functions, we can follow the same procedure and use 
the following null field instead of x'^(w) : 


dz 1 
27ri (z — wY 


(£ 

2ttiz-w 


dx 

2'Ki 

dx 


G+{x)Vo{w) 
1 


27ri (x — wY 


^G+ix)VYw) 

(Bl) 


However this procedure does not work when q = 2 since 
this is not a null field when q = 2. Instead of following 
this approach, we present a different way to obtain oper¬ 
ators annihilating the (g, 77 = 1 ) 1/9 CFT wave functions. 


using the operators already obtained in Appendixj^ The 
resulting operators are the same that would be obtained 
directly using null fields when g > 2 but this approach 
allows us to also construct a parent Hamiltonian when 
g = 2 . 

We start by multiplying the previously obtained oper¬ 
ators on the left by d-. Since d'd£ = n}\ this leads to 
new operators 


AT 
a"«"^ 


= E 7 T-’ 


,( 1 ) 


p=l,....,-2 (^^ Zj) 


= E 


(z, - ZjY 


diTi'P + d'. 


(B2) 

(B3) 


a"« = V __ d r}^^ - V 

* hi iz^ - z,Y ' * h z. - z, *■ 


(B4) 


The operator d' in A, can be replaced by 
- ^ since A^-i annihilates 

the (g ,77 = l)i wave function : 


A*’= E 


( 1 ) 


3{¥^i) 


{Zi - ZjY 


E E 

j{^i) h(^i) 


{Zi - ZhY ^ Zi - 


^ dhU^} (B5) 


A - then acts separately on H ^ and on and the op¬ 
erators AT,..., A- give zero on By keeping only 

the terms acting on and since the (g, rj = 1 ) 1/2 states 
are projections of the {q,r] = l)i states on we get 
operators annihilating the {q,r] = 1 ) 1/2 wave functions : 


A'° = E^- 

(B 6 ) 

P-1.....9-2 ( * 3) 

(B7) 

a'«- 1 -V ^ dr}^'> 

* “ ^ (z- -zY ^ ^ 



E E 

j(/i) h(iii) 


1 q nf^ - 1 (1) 

- - dhiT-l . 


{Zi - ZhY ^ - Zj 


(B 8 ) 
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